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The Tayler instability of toroidal magnetic fields in a columnar gallium 
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The nonaxisymmetric Tayler instability of toroidal magnetic fields due to axial electric currents is studied for conducting 
incompressible fluids between two coaxial cylinders without endplates. The inner cylinder is considered as so thin that 
even the limit of R ln — > can be computed. The magnetic Prandtl number is varied over many orders of magnitudes but 
the azimuthal mode number of the perturbations is fixed to m = 1, In the linear approximation the critical magnetic field 
amplitudes and the growth rates of the instability are determined for both resting and rotating cylinders. Without rotation 
the critical Hartmann numbers do not depend on the magnetic Prandtl number but this is not true for the growth rates. 
For given product of viscosity and magnetic diffusivity the growth rates for small and large magnetic Prandtl number are 
much smaller than those for Pm = 1. For gallium under the influence of a magnetic field at the outer cylinder of 1 kG the 
resulting growth time is 5 s. The minimum electric current through a container of 10 cm diameter to excite the kink-type 
instability is 3.20 kA. For a rotating container both the critical magnetic field and the related growth times are larger than 
for the resting column. 
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1 Motivation 

A known instability of toroidal fields is the current-driven 
('kink-type') Tayler instability (TI) which is basically non- 
axisymmetric (Tayler 1957; Vandakurov 1972; Tayler 1973; 
Acheson 1978). The toroidal field becomes unstable against 
nonaxisymmetric perturbations for a sufficiently large mag- 
netic field amplitude depending on the radial profile of the 
field. A global rigid rotation of the system stabilizes the TI, 
i.e. much higher field amplitudes can be kept stable. For the 
rapidly rotating regime fl > Ha (with J?a is the Alfven fre- 
quency of the toroidal field) the stability becomes complete, 
i.e. all possible modes in incompressible fluids of uniform 
density are stable (Pitts & Tayler 1985). We shall demon- 
strate in the present paper how this instability and its sta- 
bilization by rigid rotation can experimentally be realized 
with fluid conductors like sodium or gallium. There is so 
far no empirical or observational proof of the existence of 
the TI (Maeder & Meynet 2005). 

The stability problem of a system of toroidal fields 
and differential rotation has been studied in cylindric ge- 
ometry (Michael 1954; Chandrasekhar 1961; Howard & 
Gupta 1962; Chanmugam 1979; Knobloch 1992: Dubrulle 
& Knobloch 1993; Kumar, Coleman & Kley 1994; Pessah 
& Psaltis 2005; Shalybkov 2006) but in all studies only axi- 
symmetric perturbations are considered. Here as a continu- 
ation of the papers by Rudiger et al. (2007a,b) and Rudiger 
& Schultz (2010) the attention is focused to the nonaxisym- 



Corresponding author: gruediger@aip.de 



metric perturbation modes with m = 1 for real fluids. In 
particular, the possible realizations of the instabilities as an 
experiment laboratory in progress are discussed. 

2 The Taylor-Couette geometry 

A Taylor-Couette container is considered confining a 
toroidal magnetic field of given amplitudes at the cylinders 
which rotate with the common rotation rate Q. We are here 
only interested in the hydrodynamically stable regime with 
rigid rotation. The gap between the cylinders is considered 
as so wide that also a column without any inner cylinder 
can be simulated. Formally, the inner radius is 100?}% of 
the outer radius. Extreme values of fj = 0.001 are used in 
the present paper contrary to the calculations by Rudiger et 
al. (2007a) for a medium gap of fj = 0.5. The fluid between 
the cylinders is assumed to be incompressible with uniform 
density and dissipative with the kinematic viscosity v and 
the magnetic diffusivity r\. 

The container possesses the inner and outer cylinder 
with radii R m and i? ou t and with their ratio 

Rin n . 
rtout 

The radial profiles of the magnetic background fields are 
restricted for real fluids. The solution of the stationary in- 
duction equation without shear reads 

B 
R 



B 6 = AR 



(2) 
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in the cylinder geometry. A and B are the basic parameters; 
the term AR corresponds to uniform axial currents with I = 
2 A everywhere within R < R ou t, and B/R corresponds to 
a uniform additional electric current only within R < R ln . 
In the present paper we generally put 



B = 



(3) 



(see Roberts 1956; Pitts & Tayler 1985) with the con- 
sequence that the azimuthal magnetorotational instability 
('AMRI', see Riidiger 2007a) does not appear. The behavior 
of the toroidal field is thus only due to the TI for magnetic 
fields increasing outwards. It is useful to define the quantity 

Pout 



measuring the variation of B^ across the gap. Vanishing B 
leads to fi R = l/?7- In the following we fix /is = 1/fj but 
we shall vary the magnetic Prandtl number 

Pm=-, (5) 
V 

and also the (small) values of r). 

3 Equations and numerical model 

The dimensionless MHD equations for incompressible flu- 
ids are 



Re— + Re(u • V)tt = -VP - 
at 



Au + Ha 2 rot-B x B, 



Ha = 



(7) 



Here D = Ak /i?i n (i? out — i?i n ) is used as the unit of length, 
r]/D as the unit of velocity and B ln as the unit of magnetic 
fields. Frequencies including the rotation rate i? are normal- 
ized with the inner rotation rate Q m . The Reynolds number 
Re is defined as 



Re = 

v 

and the magnetic Reynolds number is 

QD 2 
Rm = . 



V 



(9) 



Sometimes it appears as useful to work with the 'mixed' 
Reynolds number 



Rm* = \/Re • Rm 



(10) 



which is symmetric in v and -q as it is also the Hartmann 
number. For Pm = 1 it is Re = Rm = Rm*. It is also 
useful to use the Lundquist number S = VPm Ha. The 
ratio of Rm* and Ha, 

Mm = ^, (11) 
Ha 

is called the magnetic Mach number of rotation. 



Applying the usual normal mode analysis, we look for 
solutions of the linearized equations of the form 

F = F(R) exp[i(kz + m<j> + ut)] . (12) 

Using Eq. (TTZl ). linearizing the Eq. (O and representing the 
result as a system of first order equations, one finds 

du R ur m 
~dR + ll +1 R U + + lku ' =0 > 



dP 
dR 



R 



ikZ+ \k + — 1 (//,• 

+ iRe(w + mf2)u R — 2J?Reit0 
- iHa 2 mAb R + 2Ha 2 Abcb = 0, 



(4) d(j) v 



dR 



R 2 



it0 — iRe(w + mJ?)it0 + 



2i—u R - 2ReQu R 



2Ha 2 Ab R + iHa^mAbrf, - i—P = 0, 



R 



dZ Z / 2 tti 
dR + R~( k + R? 



dbR 
dR 



R 



.m, 
l R k 



z — iRe(w + mQ)u z 
- ikP + iHa 2 mAb z = 0, 
ikb z = 0, 



db z i /,o m 2 \ , „ „ 1, 



1 m 



1 



OB 

Rm— = AB + Rmrotfu x B), (6) 
at 

with div u = div B = and with the Hartmann number 

B in D 



dR 



-1- Z'ft'PB - -j^mAuR = 0, 

k 2 + — J 6 - iPmRe(w + m^)6 + 
2m 

+ i-j^b R - Ru R + imAu,p = 0, (13) 



B = ^+ b 4, d4) 



where <f> u , Z and <p R are defined as 

_du£ _ du z 

0U ~ dR R ' ~ dR ' T " dR R 
and A = 1 / R- m (R ln in units of D). 

An appropriate set of ten boundary conditions is needed 
to solve the system ( fT3l >. For the velocity the boundary con- 
ditions are always no-slip while the cylinders are assumed 
to be perfect conductors. These boundary conditions are 
applied at both i?j n and i? ou t- The wave number is varied 
as long as for given Hartmann number the Reynolds num- 
ber takes its minimum. The procedure is described in detail 
by Shalybkov, Riidiger & Schultz (2002). One immediately 
finds that the sign of the real wave number k is free so that 
with the solution for k also another one with — k exists. For 
insulating walls the boundary conditions are more compli- 
cated (see Riidiger et al. 2007b). 

Tayler (1973) found the necessary and sufficient condi- 
tion 

^ (R B l) < d5) 

for stability of an ideal nonrotating fluid against nonaxisym- 
metric disturbances. Our field profile Ba — AR is thus sta- 
ble against m = and unstable for sufficiently large field 
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Fig. 1 Top: the critical Hartmann number for very wide 
container gaps (fj — > 0) between conducting and resting 
cylinders. The dotted line represents the profile dT&b with 
H = 25. Bottom: note the high quality of the representation 
( [ToT l. Valid for all magnetic Prandtl numbers. 



amplitudes, i.e. for Ha > Ha cr j t . The same is true for the 
nearly uniform fields with ps — 1- A general criterion for 
rotating flows does not exist. 

The system ( fT3l ) has been used to demonstrate that for 
resting cylinders the Ha cr it do not depend on the magnetic 
Prandtl number Pm (Riidiger & Schultz 2010). Shalybkov 
(2006) has given a similar result for m = 0. 

4 Resting container 

In the present paper the radius of the inner cylinders is as- 
sumed as very small or even vanishing (fj <C 1). We have re- 
cently shown with similar models that indeed wide gaps are 
much more suitable for TI experiments with liquid metals 
like sodium or gallium than narrow gaps. The nonaxisym- 
metric modes which are excited by the Tayler instability in 
resting columnar containers do not rotate as their drift fre- 
quencies 5R(cj) vanish. 

4.1 Critical magnetic field amplitudes 

The resulting Hartmann numbers Ha cr it for various fj are 
given in Fig. [T] They vary over many orders of magnitude 
and become very small for wide gaps (fj — > 0). For small fj 
the critical Hartmann number vanishes like 



Our calculations with the smallest R ln concern fj = 0.001 
and provide Ha cr j t = 0.00075 which leads to H ~ 25 (see 
Fig-El bottom). For given H the relation ( fT6l ) simply yields 

U (17) 



B, 



-y/pfj^vrj. 



for the critical magnetic field at the outer cylinder. Hence, 
the technical constraints to realize the TI in the labora- 
tory working with liquid metals are the following. For a 
nonrotating container with Ri n — > and i? ut = 5 cm and 
for the gallium-tin alloy used in the experiment PROMISE 
(^/pfioi/f] = 25.7 in c.g.s.), Eq. ( TFTT i leads to B ou t ~ 128 G. 
Now let I be the axial current inside the container, i.e. 

-Bout = 7-5 — ' 

O-ftout 

where R, B and I are measured in cm, G, and A. Hence, 
for i?i„ — > the total current through the fluid conductor 
becomes 3.20 kA which is surprisingly small. We have thus 
shown that it should be possible to realize the nonaxisym- 
metric current-driven TI in the MHD laboratory also with 
conducting fluids of small magnetic Prandtl number such as 
sodium and gallium. 

4.2 Growth rates 

The other important parameter for the experiment is the 
growth rate of the instability. Contrary to the critical Hart- 
mann numbers the growth rates for supercritical Hartmann 
numbers do depend on the magnetic Prandtl number. As 
it is not a mild dependence over many orders of magnitude 
of Pm a rather compact formulation is needed. We prefer 
a normalization of frequencies with the averaged diffusiv- 
ity 77* = ^/Tnj similar to the normalization of the magnetic 
field by the Hartmann number. Figure [2] shows the normal- 
ized growth rates 

" = — - < 19 > 

rj* 

for various Ha and Pm. Note that even in this formula- 
tion the growth rates strongly vary over the range Pm = 
1 . . . 10~ 6 . Fluids with equal values of v and 77 are most un- 
stable for fixed value of 77* (77* = 2.9 cm 2 /s for Gallistan, 
see below). The lower plot of Fig. |2]gives the values for var- 
ious Pm but for fixed Ha. The curve takes its maximum for 
Pm = 1. Note that the growth rates for very small mag- 
netic Prandtl numbers such as Pm = 10 -6 are two orders 
of magnitudes smaller than they are for Pm = 1. 

On the first view one takes from Fig. [2] that uj* ~ 7Ha 
which yields the relation w gr ~ 7v^^A f° r the physical 
growth rate. Here the Alfven frequency J?a is defined by 

Bout 



(20) 



Ha r 



Hi) 



1.5 



(16) 



and it is 7 < 4 for Pm = 1. Hence, for B ou t — 1 kGauss 
and an outer radius of (say) i? ut = 5 cm the growth rate is 
20 s _1 so that the growth time is 0.05 s. For Pm = 10~ 6 
(gallium) the value for 7 is only 0.05 hence the growth time 
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Fig. 2 Top: the growth rates of the instability for the rest- 
ing container normalized with the diffusion time formed 
with the averaged diffusivity 77* = *Jvrj. Fluids (or simula- 
tions) with fixed 77* have maximum growth rates if Pm = 1. 
On the right vertical axis the curves are marked with their 
value of Pm. Bottom: the same but for the fixed Hartmann 
number Ha = 4. The lo* for Pm = 1(T 6 is 0.2. It is 
77 = 0.05, hb =20. 

for i?out = 1 kG (!) with 4 s surprisingly long. Here and 
in the following the material quantities p = 6.4 g/cm 3 and 
77 = 2428 cm 2 /s for the gallium alloy GalnSn ("Gallistan") 
are used. 

However, the profiles in Fig. [2] are not linear. One finds 

that 



7^ (Ha - Ha crit ) 



(21) 



gives a better representation where 7 W depends on Pm 
rather than the value of Ha cr it- Neglecting the small values 
of the latter, 



B. 



out 



(22) 
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Fig. 3 The parameter T u in Eq. ( 1221 in its dependence 
on the inner radius 77 and the magnetic Prandtl number 
Pm (marking numbers). The Pm-dependence disappears 
for small Pm. The values hold for all containers with small 
inner radius. For Pm = 1 it is T u = 0.0009. 



is obtained for the physical growth rate. Here, 



(23) 



does no longer depend on the special (small) value of the 
geometry factor 77. Even the limit 77 — > is included. 

A weak dependence on the magnetic Prandtl number re- 
mains. Figure [3] gives the numerical values of this quantity 
for various Pm. We find the limit Pm — > as possible. For 
fixed 77 we now find fluids with Pm = 1 as most stable 
while small or large viscosities are destabilizing. 

For Pm = 10~ 6 we take from Fig. [3] the value T u — 
0.041. Hence, 



w gr = 2.0 x 10 
which yields 
500 



-3 R 2 



R2 



(24) 



(25) 



(with £>ioo = -B out /100G) as the growth time in seconds 
for the TI in gallium experiments. The size of the container 
does no longer influence the growth rate. The growth time 
for £?out = 1 kG now results to 5 s - and to ~ 2 min for, say, 
200 G. 

5 Rotating container 

Figure |4] gives for the wide-gap container the rotational 
quenching of the TI for various values of the magnetic 
Prandtl number. The lines for marginal instability in the 
Rm*-Ha cr it plane are rather straight lines. The line for 
Pm = 1 gives the maximum stabilization by rigid rotation. 
In this representation a stronger stabilization does not exist. 
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Fig. 4 The suppression of TI by uniform rotation for fj = 
0.05. The curves are marked with their magnetic Prandtl 
numbers. In this representation fluids with Pm = 1 undergo 
the strongest stabilization by rotation. The critical Hartmann 
number for resting container is 0.31 for all Pm. 

If Rm* < 4 • Ha then all fluids are unstable even under the 
influence of rotation. Note that the rotational stabilization is 
much weaker for Pm ^ 1. It is particularly weak for very 
small Pm. 

We find from Fig. |4] that for fast rotation the magnetic 
Mach number (Eq.fTTIi is almost constant for given Pm 
(with minimum for Pm = 1) but still depends on the value 
of f). We then also find that the quantity 

M = yffi Mm (26) 

does no longer depend on fj but that it still slightly depends 
on Pm (Fig. [5]). From Eq. ( f26b it results that 

fl = Mfl A (27) 

is the critical rotation rate which stabilizes the TI. Here, the 
Alfven frequency (Eq. |20l of the outer magnetic field am- 
plitude has been used. 

The question is whether the rotational stabilization of 
the TI can be probed in the laboratory. The rigidly rotat- 
ing wide-gap container with f) = 0.05 is thus considered 
also for the very small magnetic Prandtl numbers of fluid 
metals and for Reynolds numbers of order 10 3 where after 
Fig.|4]the bifurcation lines are no longer straight lines. The 
magnetic Prandtl number of sodium is about Pm = 10 -5 , 
and for gallium it is about Pm = 10~ 6 . Without rotation 
the critical Hartmann number is 0.31 for this container (in- 
dependent of Pm). With rotation the numerical results are 
given in Fig. [6] We find the rotational stabilization also ex- 
isting for conducting fluids with their very small Pm. For 
not too fast rotation the differences of the resulting criti- 
cal Ha are very small for both the fluid conductors so that 
experiments with gallium or sodium are possible. For the 
small Reynolds number of order 1000 the marginal-stable 
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Pm 

Fig. 5 The same as in Fig. [4] but valid for f) 0. Note 
that fluids with Pm = 1 can be stabilized with the smallest 
rotation rates. 
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Fig. 6 The suppression of the TI in a wide-gap container 
by rigid rotation for Pm = 10~ 5 and Pm = 10~ 6 . The 
standard Reynolds number © is given for experimental ap- 
plications. The numerical values of the outer field strength 
and the critical rotation rates are given for experiments with 
Gallistan. fj = 0.05, fi B = 20. 

magnetic field is about two times higher than for Rc = 0. 
With the density (p = 6.4 g/cm 3 ), viscosity (v = 3.4xl0~ 3 
cm 2 /s) and magnetic diffusivity (?y = 2428 cm 2 /s) for Gal- 
listan the corresponding magnetic field amplitudes (at the 
outer cylinder) in Gauss and the rotation frequencies (in Hz) 
are also given in Fig. [6] The resulting values are moderate. It 
should thus not be too complicated to find the basic effect of 
the rotational suppression of TI also in the MHD laboratory. 

Our linear code is also able to work with Pm = 0. The 
curves in Fig. [6] do not become more steep for Pm — > 0. 
There is no visible difference between the curves for Pm = 
10~ 6 and Pm = 0. 

Figure [6] does not concern to a gallium column without 
inner cylinder. However, as demonstrated above the transi- 
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Fig. 7 The parameter To for Pm = 10 5 (long line) and 
Pm = 10~ 6 (short line). 

tion from f\ = 0.05 to f\ = is easy. We take from Fig.|6]that 
close to the critical Hartmann number for resting container 

Rc = 7o (Ha-Ha crit ) 2 - 5 . (28) 

After the usual transformation to B out one finds 



with T = f? 11/4 



5/2. 



R, 



1/2 



Pm 1/4 ry 3 / 2 



(29) 



7 W . The result (1291 also holds for r) — » 0. 
Figure [7] shows the numerical values of To for the fluids 
with small magnetic Prandtl numbers. 

Using the numerical result Tq ~ 4 from Fig. |7]Eq. (|29| > 
yields 

/2~1.36B? 5 [s^ 1 ] (30) 

for a container filled with Gallistan and with i? D ut — 5 cm 
or, what is the same, / = 0.22 Sfgo as the critical rota- 
tion frequency which marginally stabilizes the TI. Hence, 
for (say) 200 Gauss the limit frequency is 1.2 Hz. which is 
much faster in comparison to the container with fj = 0.05. 

The influence of the rotation on the growth rates of the 
instability is only small. Figure [8] displays for Pm = 10~ 6 
the normalized growth rates for Re = (solid line) and 
Re = 1000 (dashed line). The stabilizing action of the basic 
rotation is obvious. 

6 The magnetic field pattern 

Without rotation the primary current-driven instability leads 
to a nondrifting steady-state solution. Isosurfaces of the pat- 
tern resulting from a 3D nonlinear simulation are shown in 
Fig.|9]for all the three field components. More details of the 
code which solves the incompressible MHD equations are 
given by Gellert et al. (2008). The boundary conditions at 
the walls are those for pseudovacuum. The container with 
fj = 0.05 is periodic in z-direction with a domain height of 
27r(i? out — i?in)- The Prandtl number Pm = 0.02 is the min- 
imum value with which the code can operate. The Hartmann 
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Fig. 8 The normalized growth rates for Pm = 10~ 6 un- 
der the influence of rigid rotation with Re = 1000 (dashed). 
The solid line gives the growth rates for the resting container 
(see Fig. |2}. Note the stabilizing action of global rotation. 

fj = 0.05, i±b = 20. 



number (Ha = 2) is supercritical with almost two times the 
critical value. The dominating modes with m = ±1 lead 
to a kidney-shaped pattern. Note the different azimuthal an- 
gles of the pattern of the azimuthal component compared to 
the others. 

From linear calculations we know that for smaller Pm 
the axial wave numbers decrease and the pattern extensions 
in axial direction become larger. 

The secondary bifurcation leads to axial oscillations 
around half the wavelength. To what extent this pic- 
ture holds for Pm = 10~ 6 is open. Because of the Pm- 
independent critical Ha and the only slightly larger Hart- 
mann number used in the simulations shown here, it is pos- 
sible that liquid metal experiments behave nearly the same. 

7 A short review 

The interplay of toroidal magnetic fields and solid-body 
rotation for incompressible fluids of uniform density fill- 
ing the gap between the cylinders of a Taylor-Couette con- 
tainer is considered. The toroidal field is the result of an 
electric axial current of homogeneous density. It is shown 
that for zero rotation the critical magnetic field amplitudes 
for marginal Tayler instability do not depend on the mag- 
netic Prandtl number. The critical magnetic field at the in- 
ner cylinder strongly depends on the gap width. It is very 
high for small gaps and it this rather low for wide gaps. For 
small enough inner radius R m the critical Hartmann number 
(of the inner field strength) runs as -R^ 5 . Note, however, that 
the critical magnetic field at the outer cylinder does hardly 
vary with the gap width (see Eq. [17). The resulting elec- 
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Fig. 9 (online colour at: www.an-journal.org) Isosurface plots of all field components. The outermost visible surface 
represents 25 % of the maximum value. Re = 0, Ha = 2, fj — 0.05, [Ib — 20, Pm = 0.02. 



trie currents necessary for TI with \m\ — 1 slightly exceeds 
3.2 kA if the material is the same gallium-tin alloy (GalnSn) 
as used in the experiment PROMISE. Such currents can eas- 
ily be produced in the laboratory. The main difficulty for 
experiments is given by the long growth times of the insta- 
bility for magnetic fields weaker than 1 kG. 

Also the rotational quenching of the TI is studied. Fig- 
ure|4]shows the rotational stabilization for various magnetic 
Prandtl numbers. To normalize the basic rotation a 'mixed' 
Reynolds number (Eq.fTOli is used in which - as in the Hart- 
mann number - the viscosities v and r\ are symmetric. The 
ratio of this Reynolds number and the Hartmann number 
(i.e. the magnetic Mach number) is free of any diffusivity. 
This ratio describes the rotational quenching of TI for var- 
ious Pm. It is shown that the most effective stabilization of 
TI happens for Pm = 1. It is weaker for both smaller Pm 
and higher Pm. Many numerical simulations with Pm ~ 1, 
therefore, could overlook the nonaxisymmetric instability 
of strong toroidal fields. 

Our calculations demonstrate that rigid rotation always 
stabilizes the magnetic field against the nonaxisymmetric 
TI. We have also shown that this rotational stabilization 
should be observable in the laboratory. Figure [6] provides 
the result that the critical Hartmann number in a wide-gap 
container of i) = 0.05 can be increased by a factor of two 
by a rigid rotation of only a few Hz. 

Acknowledgements. Many discussion with Frank Stefani and 
Gunter Gerbeth (both FZD) are cordially acknowledged about the 
theory of the current-driven instability and its experimental real- 
ization. 
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